Prom energy-density functionals to mean field potentials: a systematic derivation 
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In this paper we present a systematic method to solve the variational problem of the derivation 
of a self-consistent Kohn-Sham field from an arbitrary local energy functional. We illustrate this 
formalism with an application in nuclear physics and give the general mean field associated to the 
widely used Skyrme effective interaction. 
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I. INTRODUCTION 

The density functional theory (DFT) is one of the 
most powerful theories to deal with the intractable quan- 
tum many body problem for interacting systems with an 
arbitrary number of constituents [2]. The original pre- 
sentation of DFT from the Hohenberg-Kohn theorems Q 
establishes an exact mapping between the total ground 
state energy and the local particle density p(r) and has 
been essentially used in electron structure theory. More 
general functionals of different local one-body observ- 
ables can however be used0, @. In particular in nuclear 
physics, the so called Hartree-Fock formalism with effec- 
tive density dependent interactions [1, @, [1] can be seen 
as one of these extensions. 

The practical implementation of density functional 
theory to electronic or nuclear structure problems is done 
using self-consistent mean-field Kohn-Sham equations [9|. 
Within this framework, the interacting problem is re- 
duced to a tractable problem of non-interacting particles 
moving in an effective potential. The functional form of 
this effective potential cannot be extracted from the bare 
interactions, and different approximations are done. A 
simple local-density approximation (LDA) gives in gen- 
eral poor results, and more involved functionals including 
gradient terms, kinetic energy densities and spin currents 
are currently emploved[lCl I X II] - 

For a given arbitrary functional the definition of the 
effective mean field potential is formally straightforward 
from the variational equation 5E = tr However 

the actual computation of W is tedious and complicated 
as soon as approximations going beyond LDA are con- 
sidered. In this article, we present a systematic way to 
derive the mean field equation or equivalently the self- 
consistent Kohn-Sham equations when a general func- 
tional of any local one-body observable is considered. 



II. MEAN-FIELD FROM DENSITY 
FUNCTIONAL 

Let us consider the most general local energy-density, 
as a functional of all possible local observables such as 
densities and currents. 



We first introduce the one body density as the obser- 
vation associated to the one body density operator 



Pji 



--< > 



(1) 



where cf (cj) creates (annihilates) a particle in the single 
particle orbital \(pi) . The formalism below can be easily 
extended to treat paring by introducing the abnormal 
density Kji =< cfSj > . 

All local observables can be expressed as a function of 
the one body density as local densities p (A (r) defined 
by the relation : 



3) (r) = tr (8(y - v)ApB 



(2) 



where f is the position operator, A and B are one-body 
operators characterizing the considered density, and the 
trace, tr, runs over all single particle states. In particu- 
lar the choice A = B = 1 leads to the definition of the 
particle density, while A = B = P gives a kinetic en- 
ergy density. More examples are worked out in the next 
section. 

According to the most general version of DFT, the 
local energy density : 



E=(ij\H\ifj) = J H(r)dr 



(3) 



can be expressed as a functional of an ensemble of local 
densities {p A)B (v)} associated to an ensemble of opera- 
tors pairs (A, B) such as: 



(4) 



Using the variational derivation of mean field theory Q , 
the mean field hamiltonian is defined by the equation 

5E = tr (Wdpj (5) 

We can thus write the differential of the total energy in 



2 



the form 



SE = J 8H{r)&r 



(A,B) J {A ,B) 

on 



E 



111 E a^~^ A5pB 



(r)tr ((5(f - r)A5pi?) dr 



(6) 



Identifying expression ([6]) with the mean-field hamilto- 
nian definition ([5]) we get : 

w= E ^,b)= E ^^r^ 1 ^ 

(A,B) (A,B) UtJ iA ,B) 



A. Examples of usual operators 

In this part, we will identify the various terms, den- 
sities and currents, which are used in the most common 
density dependent energy functionals with specific local 
densities p (A B) (r) 

1. Particle densities 

Let's start with the total particle density p(r) that we 
can write at first : 



p( r ) = E( r ' a i^i r ' a ) 



(8) 



where the sum over a runs on all additional quantum 
numbers such as spin and isospin. Introducing the oper- 
ator 5 (r — f ), we can rewrite this expression in a trace 
form: 

p(r) = tr (6 (r - f ) p) (9) 

showing that A = B = 1 i.e. 

p(r)= P(11) (r) (10) 

Nuclear functionals generally use spin and isospin den- 
sities 



Pta(r) =P (n ,„ a) (r) 



(11) 



where ft is the t component of the isospin Pauli operator 
and (7 S the s component of the spin Pauli matrix. Of 
particular importance for asymmetric nuclear matter and 
exotic nuclei is the isovector particle density 



where T3 is the third component of the isospin operator. 

Extending the f and a operators to include the iden- 
tity To = 1, and oo = 1 the definition includes the 
particle density (fT0|) and the isospin density (| 1 2[) . and 
can be taken as the most general. Note that since spin 
and isospin commute the order of the operators T and 
a in pip is not important, pt s — pst- It should also be 
noticed that, in the nuclear case, the functional is often 
written in terms of proton and neutron density. In such 
a case one should use the projector II„ — (1 + t^)/2 for 
neutrons and Il p = (1 — f^)/2 for protons. 



2. Kinetic densities 

The case of the kinetic energy density T(r) can be 
treated in the same way : 

nr) = EE< r > a i pJp.m (1 3 ) 

a k 

showing that T(r) is the sum of 3 densities with 



T (r)= E ^» 



(14) 



fe=l,3 



This definition can be extended to introduce spin and 
isospin as in (jlip . in particular the isovector kinetic en- 
ergy density 



T 3 (r) = V p. ,(r) 



(15) 



fe=l,3 



More generally the spin and isospin kinetic energy den- 
sities can be defined by 

TUr) = P lPk , w „Jr) (16) 



3. Spin currents 

Another important quantity in density functional the- 
ory is the spin current J(r) which is needed to describe 
the spin-orbit interaction. It is defined by : 

J(r) = -i ^2 [V r <p n (r,a)]ipn(r,a') x {a'\a\a) + c.c. 

n,a,et' 

where a is the Pauli matrix. Such term can be easily 
derived from the spin and isospin currents 



Jkts(r) = p at) {r) + c.c. 



(18) 



fls(r) = % 3ll) (r) 



(12) 



using J,(r) = J2ks=i e ^ JkOs(r), where e is the anti- 
symmetric Levi-Cevita tensor. 
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4- Gradient terms 



with 



Energy functional arc often written in terms of gradi- 
ent of densities p (A B) (r) . Such terms can also be viewed 
as (combination of) generalized densities p , B/) (r) using 
the relation 



* v P(A.B)( r ) = tr^(r-r) 



p,A P B 



.Bp) 

(r) (19) 



Higher order gradients can be similarly expressed if the 
above expressions are recursively applied: 



3=0 



tr 



(8(r - r)p j ApBp n - 



(20) 



3=0 



where C„ are binomial factors. 



a i 

a-> 
a 6 
a 7 
as 
a<:> 
an 



^(2 + xo), a 2 = -| (2x + 1) 
^(2 + z 3 ), 04 = ~ (2S3 + 1) 
g [*i (2 + ar 1 )+i a (2 + a*)] 
~ [i 2 (2x 2 + l)-*i (2a:i + 1)] 



— [3ti (2 + ari)-t2(2 + a a )] 

-^[*2 (2x 2 + l) + 3ti (2£i + l)] 

Wo 1 , 

' a w = -r^{tiXi+t 2 X2j 
2 lo 

1 



16 



(h - t 2 ) 



We first write the various terms in terms of generalized 
densities p (A B) 



K 
Ho 
H 3 



v^3 

2m 

aipl A) +a 2 (p( nnjl) +P(n pil) ) 

2+ a | a / 2 , 2 \ 

°3P(iTi) + a 4P(i,i>cn„,i) +P(n p .i)) 



III. APPLICATION TO SKYRME 
FUNCTIONALS 



As a practical application of the above formalism, let 
us consider the Skyrme nuclear energy-density functional 
" G3 defined by : 



We// = a 5EP( Pfc , Pfc ,P(i,i 



fe=l 

a6 ^(n„p fc ,p fc )P(n„,i) + P(n pPfc , Pfc )P(n p ,i)/ 
3 

2 



ft/in = -^Xm^,!) -p ihPh) y 



k=l 



) f(n„,j>,,v 



W = /C + Ho + 7^3 + He// + W/i n + Hso + (21) 

where the coefficients are linear combination of the 
usual Skyrme parameters: 



n so = ia g £iks[p (Pk>trs y(p (p . A) -P( llPi) ) 

iks — l 

+ P(n pPk ,<T S ) ■\P(,pi,n p ) ~ P(n p , Pi )) 
^<n„ Pfc , CTs ) -(P( Pi .n„) — P(n„, Pi ))i 



K 


T 
2m 




Ho 


= aip 2 + a 2 (pp + p 2 n ) 




H 3 


= p a [a a p 2 + a 4 (Pp + pl)\ 






= a 5 Tp + a 6 [T p pp + T n p„ 




Hfin 


= a 7 (\7p) 2 + a 8 [{V Pp ) 2 + 


(Vp„) 2 


H so 


= a 9 [J.Vp + Jp.Vpp + J„ 


•Vp„] 


Hsg 


= amJ 2 + an [Jjj + J 2 ] 





^iks^ik's'P 



iksk' s' — 1 
3 



+ an Siks^ik' s' (P( Pfc ,n p cr s ) P( Pfc , ,n p <T s ,) 

iksk' s' — 1 

^( Pfc .n„ CTs )P( P ,.,.n„ CT ,)) 



Using relation (0), the involved functional derivatives 
of eq.© are reduced to simple derivatives of scalar ex- 
pressions: 



4 



while the central potential reads 



K 2m 

W Ha = 2aip(f) + 2a 2 (fl n p rl (f) +fl pPp (r)) 

Wn 3 = (2 + a)a 3 p 1+a (r) 

+ aa^-\v){pl{v)+pl{v)) 

+ 2a i p a - l {v){Ii nPn {v) + n p p p (f)) 



fe=i 

3 

+ ^2 a 6 fl n p k p n (r)p k + a 6 TL n T n (r) 
fe=i 

3 

+ X! a 6n p p fe p p (f )p k + a e fl p T p (r) 



k=l 



n fin 



-2a 7 Ap(f) 



- 2a 8 (n„Ap n (f) + n p App(f )) 
W Ws< , = a 9 (Vp(f) + VJ(f) 

+ n„Vp„(f) + n„vj„(f) 
+ n p Vp p (f) + n p vj p (f)) 

W Hsg = 2a 10 J(?) 

+ 2on(n n j n (f)+n p jp(f)) 

The one body mean field operator for the charge q = 
n,p can be expressed in a more compact form as : 



Wq = Pl 2Mff + Uq{i) + mq{ * ] + Yqim x 5 (22) 

where the effective mass is 



2m(f) 



— +a 5 p(r) + a 6Pq (r) 



(23) 



U q (r) = 2a lP (r) + 2a 2 p q (r)+a 3 (2 + a)p 1+a (r) 

+ a i ap a -\v){pl{v)+pl{v)) + 2a i p q {v)p a {v){2A) 
+ a 5 T(r) + a 6 T q (r ) 

- 2a 7 Ap(f) - 2a 8 Ap,j(f) + a 9 (vJ(f) + VJ,(?)) 
and the spin orbit part is given by 

V,(r) = a 9 (Vp(f) + Vpq(f)) + 2ai J(f) + 2onJ,(f) 

(25) 

consistently with the expressions given in the litera- 
ture. 



IV. CONCLUSION 

In this paper we have introduced a simple and pow- 
erful method to solve the variational problem of the 
derivation of a mean field potential from and arbi- 
trary energy-density functional. Observing that the 
most general DFT energy functional can be expressed in 
terms of generalized local one-body densities p (A B) (r) = 

tr (<5(f - r)ApB) , the functional derivative defining the 

Kohn-Sham or HF field W = SE/8p can be reduced to 
the calculation of ordinary derivatives dTi/dpr^ B) of the 
associated energy density. As an example we have worked 
out the well known expression of the general Skyrmc 
mean-field. This method gives a simple technique to take 
into account new terms in future more sophisticated gen- 
eralized functionals. The autors want to thank M. Ben- 
der for the useful discussions. 
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